Universality in bosonic dimer-dimer scattering 
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Bosonic dimer-dimer scattering is studied near the unitary limit using momentum-space equations 
for the four-particle transition operators. The impact of the Efimov effect on the dimer-dimer scat- 
tering observables is explored and a number of universal relations is established with high accuracy. 
The rate for the creation of Efimov trimers via dimer-dimer collisions is calculated. 

PACS numbers: 34.50.Cx, 31.15.ac 



I. INTRODUCTION 

Few-body systems with large two-particle scattering 
length a possess a number of universal properties that are 
independent of the short-range interaction details. The 
three-particle system is simple enough such that analytic 
or semi-analytic derivation (at least of some) of those 
properties is possible [l| . Probably the best-known exam- 
ple is the Efimov effect, the existence of an infinite num- 
ber of weakly bound trimers with geometric spectrum in 
the unitary limit o — >• oo [2j. The four-particle system, 
however, is too complicated for analytic approaches and 
therefore so far has been investigated mostly numerically. 
While the bound state properties have been calculated 
by several methods 043) the description of the four- 
particle continuum still constitutes a serious technical 
challenge, especially in the universal regime with many 
open channels; first studies have been made using the 
hyperspherical formalism [8], (9| ■ Given the success of the 
four-nucleon scattering calculations [10|, [HI based on ex- 
act Alt, Grassberger, and Sandhas (AGS) equations [12| 
that were solved in momentum-space, we applied them 
recently also to bosonic atom-trimer scattering [13, [l4( 
and established a number of universal relations between 
observables. In the present work we extend the technique 
of Ref. 13] to study the universal physics in the bosonic 
dimer-dimer scattering; in our model only one (shallow) 
dimer state exists. We note that this problem has al- 
ready been considered in Ref. [15| using the adiabatic hy- 
perspherical framework. However, the AGS method en- 
ables us to determine the universal properties with much 
higher accuracy; in particular cases there are significant 
corrections to the results of Ref. [15J . Furthermore, we 
also present more detailed analysis of dimer-dimer scat- 
tering observables. 



II. DIMER-DIMER SCATTERING EQUATIONS 

A rigorous description of the four-particle scatte ring 
can be given by the Faddeev-Yakubovsky equations [16| 
for the wave-function components or by the equivalent 
AGS integral equations [l2J for the transition operators. 
The latter ones are better suitable for the multichan- 
nel scattering problem where the on-shell momenta and 



binding energies of the asymptotic states differ by many 
orders of magnitude |17| . Furthermore, the AGS equa- 
tions lead more directly to the observables since the on- 
shell matrix elements of the transition operators calcu- 
lated between the components of the corresponding ini- 
tial and final channel states yield scattering amplitudes 
18]. For the system of four identical bosons we use the 
symmetrized form of the AGS equations. In the case of 
the atom-trimer scattering they are given in Ref. [13], 
whereas for the dimer-dimer scattering they read 

^12 = (GotGo) + Pz4:UiGatGaUi2 + C^GoiGoZ^, 

(la) 

U 22 = (1 + P 34 )U 1 G Q tG U 12 . (lb) 

The transition operator U22 describes the elastic scatter- 
ing and U12 the transfer reactions leading to the final 
atom-trimer states. The dynamic input is the two-boson 
potential v. It yields the two-boson transition matrix 



t = v + vGot 
and the symmetrized transition operators 
U a = P a G + P a tGoU a 



(2) 



(3) 



with a = 1 and 2 for the 1 + 3 and 2 + 2 subsystems, re- 
spectively. In the above equations Go = (E + iO — Ho) -1 
is the free resolvent for the four-boson system with en- 
ergy E and kinetic energy operator Hq. The permutation 
operators Pij of particles i and j and their combinations 
Pi = P12 P23 + P13 P23 and P 2 = P13 P 2 4 are sufficient 
to ensure the correct symmetry of the four-boson system 
when a special choice of the basis states is used [13> LLSj ■ 
Namely, the basis states must be symmetric under ex- 
change of particles 1 and 2; in addition, for the 2+2 sub- 
system they must be symmetric also under exchange of 
particles 3 and 4. 

We employ momentum-space partial-wave framework 
[13 ] to solve the AGS equations (P). After the dis- 
cretization of momentum variables the system of integral 
equations becomes a system of linear algebraic equations 
whose dimension is significantly reduced by using a sepa- 
rable two-boson potential v = \g)X{g\ of rank 1 acting in 
the S wave only. Nevertheless, higher angular momen- 
tum states with l y ,l z < 3 (see Ref. pJl for the notation) 
are taken into account for the relative motion in 1 + 2, 



A Energy 




FIG. 1. Schematic representation of the four-boson en- 
ergy spectrum; only one family (n) of multimers is shown. 
The tetramer (T) energies, atom-trimer (a+t), dimer-dimer 
(d+d), and dimer-atom-atom (d+a+a) thresholds are dis- 
played as solid, dashed, dashed-dotted, and dotted lines, re- 
spectively. They intersect at special values of a described in 
the text. For a better visualization only qualitative relations 
are preserved. 



1 + 3, and 2 + 2 subsystems such that the results are con- 
verged to at least four-digit accuracy with respect to the 
partial-wave expansion. 

To prove that our results are indeed independent of 
the short-range details of the interaction, we use two 
different choices of the potential form factor (k\g) = 
[l+c 2 (fc/A) 2 ]e-( fe / A ) 2 with c 2 = and c 2 = -9.17. It was 
shown in Ref. |13| that they yield a very different trimer 
ground state and related observables, but, nevertheless, 
lead to the same universal relations for the atom colli- 
sions with highly excited trimers where the finite-range 
corrections become negligible. 



III. RESULTS 

We study the dimer-dimer scattering observables as 
functions of the atom-atom scattering length a. We fix 
the geometry of the potential, i.e., the cutoff parameter 
A, whereas the strength A is adjusted to reproduce the 
desired value of a. We present the results as dimension- 
less ratios that are independent of A and particle mass 
rn. 

A schematic a-dependence of the energy levels and 
thresholds in the four-boson system is shown in Fig. Q] 
As a reference point for a we choose the intersection of 
the dimer-dimer and nth atom-trimer thresholds, i.e., 
b n = 2bd at a = a n d , where bd and b n are the binding 
energies of the dimer and the nth trimer, respectively. 
A point important for three-boson physics is a = a^ 
where the nth trimer is at the atom-dimer threshold, i.e., 
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dd 1 dd 


dd I dd 

a n,2l a n 


at d /a d n 


atV mb n 


3 

4 


0.32335 
0.32352 


0.99948 
0.99947 


6.8019 
6.7896 


0.070601 
0.070754 


3 

4 


0.32369 
0.32354 


0.99948 
0.99947 


6.7809 
6.7886 


0.070857 
0.070768 



TABLE I. Special values of atom-atom scattering length cal- 
culated using potential form factor with C2 = (top) and 
c 2 = -9.17 (bottom). 



b n = bd- Furthermore, the existence of two tetramers 
(k = 1,2) associated with each Efimov trimer was pre- 
dicted in Refs. [a, 0]- Except for the lowest two, all other 
tetramers lie above the trimer ground state (n = 0) and 
therefore are unstable bound states [13, LL4| with finite 
width r„ j.; their position relative to the four free par- 



ticle threshold is — B„ &. At a = 



- dd 
l n,k 



the correspond- 



ing tetramer intersects the dimer-dimer threshold, i.e., 
B n ,k ~ 2&d, Tn.fc = 0, leading to resonant effects in 
the low-energy dimer-dimer scattering. The relations be- 
tween these special a values are collected in Table Q] for 
n = 3 and 4; to achieve the universal limit with good 
accuracy we concentrate on the regime with n and a 
large enough, aA/2 ranging from 10 3 to 10 6 . The dif- 
ferences between the results obtained with c 2 = and 
c-2 = —9.17 in the potential form factor are below 0.3% 
in the n = 3 case and below 0.015% in the n — A case. 
In addition we include a n y/mb][ where b^ is the nth 
trimer binding energy in the unitary limit. The agree- 
ment between our predictions and the semi-analytical re- 
sult lim^oo af lV / mb~^ « 0.0707645 given in Ref. Q is of 
the same quality as for other ratios in Table HI i.e., it is 
better than 0.015% in the n = 4 case. We therefore con- 
clude that our n = 4 calculations provide the universal 
results for dimer-dimer collisions with high accuracy, 



ai d Ja d n d = 0.3235(1), 

< 2 /«f - 0.99947(2), 

a d n d /ai = 6.789(1); 



(4a) 
(4b) 
(4c) 



the errors are estimated from the residual dependence 
on c 2 and included l v , l z . The above ratios calculated in 
Ref. [H with local potential at n = 2 are 0.352, 0.981, 
and 6.73, respectively. Thus, they differ from our results 
by 1 to 9%. In particular, our l—a^ 2 / a rf value is smaller 
by a factor of 35, indicating that the shallow tetramer is 
much closer to the corresponding atom-trimer threshold 
and thereby much stronger affects the atom-trimer low- 
energy scattering jl4| . 

In the following we do not demonstrate explicitly the 
convergence of our results, however, they are checked to 
be independent of c 2 and n for n > 3 with good accu- 
racy. In Figs. [2] and [3] we show the dimer-dimer scatter- 
ing length Add and the effective range parameter r^ as 
functions of a. In the universal limit Addja and Tdd/a 
are log-periodic functions of a, i.e., the behavior shown 
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FIG. 2. (Color online) Real and imaginary parts of the dimer- 
dimer scattering length as functions of the atom-atom scat- 
tering length. In the left inset KeAdd/a (solid curve) and 
ImAdd/a (dotted curve) are shown in the vicinity of the (n, 1) 
tetramer using linear scale. In the right inset ImAdd/a, is 
shown in the vicinity of the dimer-dimer and atom-trimer 
threshold intersection. 



in Figs. [5] and [3] is repeated when a increases by the Efi- 
mov factor e" /s ° « 22.694 [1]. Since A dd /a and r dd /a 
vary over many orders of magnitude, we use logarithmic 
scale but distinguish between positive and negative val- 
ues. A dd exhibits a typical resonant behavior at a = a^ k ; 
for k = 1 it is shown with finer resolution in the left in- 
set of Fig. [2j At the intersection of the dimer-dimer and 
atom-trimer thresholds A dd has a cusp that is especially 
pronounced for the imaginary part as demonstrated in 
the right inset of Fig. [5] 

The validity region of the effective range expansion 
is considerably narrower than the difference between 
the thresholds \2b d — b n \. Thus, the dimer-dimer ef- 
fective range parameter r dd is ill-defined at the inter- 
section of the dimer-dimer and atom-trimer thresholds; 
Rer^d (Imr dd ) diverges when a approaches aff from be- 
low (above) as demonstrated in Fig. [3] In contrast to 
A dd , r dd shows no resonant feature at a = a d f k but ex- 
hibits very rapid variations when KeA dd goes through 
zero, i.e., at a/aff w 0.200 and 0.951. In a very narrow 
region around these a values Im r dd has positive and neg- 
ative peaks whereas Re r dd has two positive (negative) 
peaks separated by a much more pronounced negative 
(positive) peak at a/aff w 0.200 (0.951); see the inset of 
Fig. [3J However, despite the dramatic variations of r dd , 
the quantity A^r^l remains almost constant in these 
regions; this is consistent with the corresponding behav- 
ior of the two-particle system [l| . 

Next we study the dimer-dimer complex phase shifts 
6l and inelasticity parameters tjl = e - 2ImA £ a s functions 



FIG. 3. (Color online) Real and imaginary parts of the dimer- 
dimer effective range parameter as functions of the atom- 
atom scattering length. In the inset Re rdd/a (dash-dotted 
curve) and Im rdd /a (dotted curve) are shown in the vicinity 
of ReAdd = using linear scale. 



of the relative dimer-dimer energy E dd up to the atom- 
atom-dimer threshold; L is the orbital angular momen- 
tum for the relative dimer-dimer motion. S- wave (L = 0) 
results depend strongly on a as shown in Fig. 3] for few 
specific a values. We remind that the standard effec- 
tive range expansion is not valid in the case of coinciding 
thresholds a/aff — 1. A very rapid decrease of Re<5o at 
a/a^ d = 1 is caused by the proximity of the (n, 2) un- 
stable tetramer. At a/aff < 1 the cusp in Re <5o and the 
descent of 770 correspond to the opening of the nth atom- 
trimer channel. Due to centrifugal barrier these features 
are absent in the D wave (L = 2). Re 82 are small but, 
depending on a, can be of both signs, while 772 deviate 
from 1 quite significantly as E dd increases. Thus, near 
E dd ~ b d the D-wave contribution to the inelastic cross 
section may exceed the one of the S wave, the ratio be- 
ing 5(1 — r/|)/(l — ?7o)- Higher waves can be neglected in 
the considered energy regime E dd < b d , i.e., 5l f» for 
L > 4. 

The cross sections for the dimer-dimer elastic scatter- 
ing a(dd -> dd) and transfer/rearrangement reactions 
a(dd —> n') leading to an atom plus trimer in the n'th 
state are shown in Fig. [5] At very low energy <r(dd — > dd) 
become nearly energy-independent (this is true also for 
a/a d n d = 1 at E dd /b d < 10~ 5 ) while a{dd -> n') in- 



crease with decreasing E dd nearly as E 



-1/2 

dd 



The 



ccption is a(dd — > n) in the case of coinciding thresholds 
a/aff = 1 where the increase slows down until cr(dd — > n) 
becomes almost constant at E dd /b d < 10~ 5 . Thus, the 
very low energy behavior of the cross sections is consis- 
tent with the Wigner threshold law. The reason why for 
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FIG. 4. (Color online) Energy-dependence of the S- and D-wave phase shifts and inelasticity parameters for the dimer-dimer 
scattering at specified values of the atom-atom scattering length. For a better comparison we set Re^L = at Edd = 0. 



a/aff = 1 the region of nearly constant a(dd — > dd) 
and a(dd — > n) is so narrow is the proximity of the 
(n, 2) unstable tetramer. The inelastic scattering, apart 
from the very narrow regions just above the atom-trimcr 
thresholds, is dominated by the open channel with the 
shallowest trimer; the reactions leading to more deeply 
bound trimers are strongly suppressed. The elastic cross 
section appears to be remarkably sensitive to the atom- 
atom scattering length. For example, a/aff = 0.814 cor- 
responds to Kerdd ~ leading to an almost constant 
a{dd -)• dd) at E dd /b d < 0.1 while for a/af = 0.951 
where Re Add ~ the elastic cross section is very small 
at Edd — but then increases with energy, in contrast 
to the other shown cases. Furthermore, the opening of 
the nth atom-trimer channel affects the elastic dimer- 
dimer cross section in opposite ways for a/aff — 0.814 
and 0.951. 

The above results are indispensable for the description 
of the dimer-trimer conversion process, i.e., the creation 
of Efimov trimers via rearrangement reaction in a gas of 
dimers. The total trimer creation rate in the dimer-dimer 
scattering j3dd = J2n' @dd nas contributions from all open 
atom-trimer channels, 



Pdd = (vdd<r(dd->n')), 



(5) 



where Vdd — y^Edd/rn is the relative dimer-dimer ve- 
locity and (...) denotes the thermal average. At van- 
ishing temperature the trimer creation rate Pdd(0) = 
— (8Tr/m)ImAdd is determined by the results of Fig. [2] 



thus, it gets resonantly enhanced around a 



,dd 



The 



results at finite temperature T, assuming the Boltzmann 
distribution for the relative dimer-dimer energy, are pre- 



sented in Fig. [5] for few selected values of a. The shown 
examples differ remarkably in their T dependence: at 
very low temperatures f3dd decreases with T very rapidly 
at the resonance (a/a„ = 0.3235) while in other cases its 
variation is considerably slower. A more pronounced in- 
crease of fidd at higher temperature for a/a^ d = 0.814 and 
0.951 is due to the opening of the nth atom-trimer chan- 
nel. We do not show predictions for individual /3% d , but 
the results of Fig. [5] clearly indicate that shallow trimers 
are produced most efficiently. 

Finally we note that f3dd was calculated also in Ref. [15[ 
where it was called the dimer-dimer relaxation rate. 
There is a good qualitative agreement between our f3dd 
results and those of Ref. [15j. There are, however, some 
quantitative differences, e.g., at a/aff — 1.5 our fidd de- 
creases with T slower. At least to some extent this is 
probably caused by the significant D-wave contribution 
taken into account in our calculations. 

The only available experimental data [20j refer to the 
regime where a exceeds the interaction range by a factor 1 
to 10 and therefore non-negligible finite-range corrections 
can be expected. In this respect our calculations present 
no improvement over those of Ref. [15] that could provide 
only qualitative explanation of some features of the data. 
However, our predictions would be valuable for the future 
experiments performed in the universal regime. 



IV. SUMMARY 

We studied bosonic dimer-dimer scattering near the 
unitary limit. It is a complicated multichannel four- 
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FIG. 5. (Color online) Elastic and transfer cross sections for the dimer-dimer scattering at specified values of the atom-atom 
scattering length. 
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particle scattering problem involving, in the present cal- 
culations, one dimer-dimer and up to five atom-trimer 
channels with very broad range of binding energies. Ex- 
act AGS equations were solved using momentum-space 
techniques. The log-periodic structure of the dimer- 
dimer scattering observables was found to be independent 
of the short-range potential details. Universal results for 
the dimer-dimer scattering length, effective range, phase 
shifts, elastic and transfer cross sections, some tetramer 
properties, and trimer creation rate were obtained with 
an accuracy considerably higher than in previous works. 



FIG. 6. (Color online) Temperature dependence of the trimer 
creation rate at specified values of the atom-atom scattering 
length; a/a n = 0.3235 corresponds to the (n, l)st resonance 
that leads to a peak of 0dd at T = 0. 
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